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ABSTRACT 


This thesis is a study of the Omega navigation system as 
currently implemented in the P-3C aircraft. The possibility 
of using a microcomputer to solve the internal processing 
functions is investigated. Data flow graphs were applied to 
the velocity and navigation processing function in the Omega 
system. These graphs assisted in the development of the 
PL/M code which implements the function. The four PL/M 
subroutines that were written can compute the velocity and 
navigation equations in sufficient time and with sufficient 
qO G umay ЕО encourage additional research InCo a 
microcomputer implementation of the remaining internal 


functions of the Omega systen. 
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This thesis examines the use of a microcomputer to solve 
the internal processing functions of the current P-3C Update 
Omega navigation system. The most important reason for 
studying &a microcomputer for the Omega navigation system is 
that during an excessive system load in the P-3C aircraft, 
the Omega signal processing 1s partially inhibited. 
Operating in this inhibited mode appeared unnecessary in a 
sophisticated system such as the P=3C incorporates, 
especially with the existing technology of microprocessors 
and distributed systems (Refs. 1, 2, апа 31. with a 
microcomputer dedicated to Omega navigation procsssing, this 
inhibited mode of operation could be eliminated and position 


updating information could be available continuously. 


The analysis of the existing velocity and navigation 
processing routine utilizes the current equations in the 
published manual on the P-3C Update system functional 
description for Omega. The main tool for the analysis is 
the data flow graphs of the existing equations. This gives 


efficient and optimal PL/M code. 


бесінен» Id» of this thesis presents fundamental 
information applicable to the worldwide Omega navigation 
system and the AN/ARN-99(V) receiver-converter. Section III 
introduces the various internal processing functions of the 
system and presents a data flow graph depicting the 
parameters which apply to each of the functions. Section IV 
presents the analysis of the velocity and navigation 
processing routine with applicable data flow graphs. 


Section V gives the conclusion of the study. 





OMEGA is a worldwide, very low frequency (VLF), radio 
navigation system which utilizes the radiation from eight 
transmitting stations to provide global coverage to 
aircraft, snips, land-vehicles, and submarines for accurate 
and reliable positioning. Table I lists the presently 
available transmitting stations with appropriate letter 
designators, coordinates, and direction cosine values [Ref. 
4]. The position of each transmitting station is shown in 


Figure 1. 
. SYST8W DRINCIPLES 


Each of the eight transmitting stations radiates 
continuous, sinusoidal wave bursts at 10.2 kHz, 13.6 kHz, 
and 11.3 kHz. These signals are phase locked and time 
synchronized to Universal Time such that ail three signals 
start at zero vaiue with positive slope at 0000 hours 
Greenwich Mean Time (GMT) and repeat at ten second intervals 
[Ref. 5]. Figure 2 depicts a standard ten second interval 


transmission pattern. 


In order to understand the complexity of the Onega 
navigation system and the operations performed by the 
receiver and the computer software functions, a simplified 
view of the systems operation is presented first. Assume 
some arbitrary transmitting station is radiating on one of 
the three frequencies, and, at Some distance d, a receiver 


Ше dad Gal 1 ring the signal for future processing. The 





STATION LOCATION LAT/LONG BERECTEON COSINE 


Norway O0 R251 T2. 39 N +99135351329 
130 8 12.65" +0.39172334 

| +0.09 153733 

Liberia N +0.98221041 
к +0.98117517 

-0.16588382 

Hawaii 16.90" N +0.36 2833231 
52.70" W -0.86296869 

-0.35162108 

57.20! +0. 721443222 


08.77" W -0.10039077 
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Island 2029! *T05 5321403595 
+0.76823587 

Argentina 12.53! -0.68022941 
27.60" Y +0.30756230 

-0.66535207 

Trinidad 06.20" +0.18449565 












+0.46687116 
-0.86486570 
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-0.52409936 
+0863633039 
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OMEGA SYSTEM TRANSMITTER LOCATIONS 
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transmitted signal appears as a standard sinusoidal wave. 
Assuming no attenuation of the signal, a perfect waveguide 
formed ov the earth and the ionosphere for the signal to 
ООО chrouga, no interfering carrier or noise signal, 
and the transmission traveling at the speed of light, the 
received signal would appear as a sinusoidal wave but offset 
at the receiver in time t by d/c, the time necessary for the 


Signal to travel distance d at the velocity of light c. 





t 1 

а/с 
O E ЙЫ) m0 O<t<d/c 
SIN (wt) t20 SIN w(t-d/c) t»d/c 


However, this ideal situation is rarely achieved. The 
propagated signal is affected by several factors described 
Ша SubSesHenct section. Due to the repetitive nature of 
sinusoidal waves, the precise number of wave lengths nust be 
computed between the transmitter and the receiver. Also, 
the absolute time t of the transmitter must be known. These 


Пас аге calculated in the combinational filter. 


Dee) fi wng With the use of OMEGA signals is 
accomplished by two methods. The first method utilizes 
phase difference information between stations and the 
resultant hyperbolic lines of position. The second method 
uses circular lines of position obtained by measuring phase 


with respect to a stable frequency standard. The P-3C Omega 
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Mayigation system uses the circular lines of position or 


rho-rho metnod to deternine position fixes. 


Of primary importance in utilizing the OMEGA VL? signals 
for the rho-rho method are the ability to detect and measure 
ре рпасе оі а received signal, and the ability to predict 
the phase of the signal. In order to detect and measure the 
phase of a received signal using the rho-rho method, the 
first process performed in the Omega navigation system is to 
er are he internal clock of the computer with the 
transmission purst pattern of the stations. This 
synchronization process is describeä in more detail in a 
subsequent section. The overall effect allows the computer 
program to know what station and frequency is being 
processed at any internal system tine t. If all three 
frequencies £rom all eight stations could be received at any 
ООПТ on the earth's surface, the ten second reception 
pattern would appear aS in Figure 3. For a more realistic 
reception pattern, Figure 4 depicts the three closest 
transmitter stations to Hawaii, and Figure 5 represents the 


signal pattern an aircraft would receive at Hawaii. 


Once synchronization is compieted, the Omega navigation 
system ccmmences to process the received signals througn 
Various scitware filters in order to measure the phase. The 
rilters which process the received signal are composed of 
three burst fllters (ie. one for each frequency) and 24 
tracking filters (ie. one for each of the tnree frequencies 
C por he eight stations). Тһе primary function cf the 
burst filters is to calculate the phase measurement and the 
phase variance for each frequency and station received by 
the receiver. The tracking £ilter, on receiving the phase 
Beasurement and variance of the phase measurement from the 
burst filter, utilizes these values to acquire a more 
accurate estimate of phase and phase variance over several 


eens second time intervals. 
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Bu etracking filter derforms two different updates on 
F Sor SU Filter data. Тһе first update is a time update. 
This update is required because of the displacement of tae 
alrcraft receiver between successive measurements of the 
received signal. The time update procedure projects the 
last estimate of phase through time so that it reflects the 
new position of the aircraft. These procedures are defined 
as rate aiding and are dependent on velocity sources 
external to tne Omega system. In estimating the change in 
position, it is possible to estimate what the cnange should 
ре Бәке а pantuicular station and frequency. This predicted 
phase change is then added to the last estimate of phase to 


create a new phase estimate. 
ли (т) = DIRE (1-1) + СГА фр (Е) + AQERR (t-1) ] = At 


where Д фра (ї) іѕ tne phase rate which represents the 
average rate of change of phase along the arc connecting a 
Eon and the aircraft, and A DERR(t-1) is the estimate of 
the error in the pnase rate and is calculated in tae 
neasurenent update Sect lon. The second update is a 
measurement update, and is performed every time the burst 
filter supplies a measurement of phase and phase variance to 
M tracking filter. This measurement update procedure 
combines the rate aiding phase estimate from the tracking 
filter with the burst filter phase measurement to produce a 
bester estimate of the phase based on tne  foliowing 


equation. 

IN ТЕКТЕ) + [ AQDR(t) +A DOERR (t) ] * (9.425) 
where ф іс defined to be the angular difference of phase. 
The combination of the phase measurement from the burst 


filter aná the phase estimate from the tracking filter is a 


4eighted average, where tne coefficients are computed using 
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Piempidee Vatiance From the burst filter (c ?0JK) and the 


phase variance rrom the tracking filter (o ?0TRK). 


ФТЕК (2) = OTRK(t-1) + ARCTAN (o 20TRK * S) 
( (с 2JK + o 20TRK) * C) 


Mete. S = SIN p andá C = COS 0. 


Once the measurezent of phase nas been calculated, it is 
supplied to the combinational filter for processing along 
Ne a phase estimate calibrated in the propagation 
predicticn procedure. With the measured phase and the 
predicted phase available for position fixing, the last 
requirement to ое implemented for utilization of the rho-rho 
technique of position fixing is to synchronize the receiver 
О-ОО атт with the transmitter oscillator. This is 
accomplished in the combinational filter by computing the 
time difference between the transmitting station phase and 


the receiver oscillator phase. 
1. Signal Coverage 


The VLF signals utilized by the OMEGA system are 
transmitted cver extremely long distances or tha earth's 
surface, because tney are propagated through a natural wave 
guide formed by the surface of the earth and the ionosphere. 
Due to the advantageous properties of high phase stability 
and low attenuation rates or VLF signals, no modulation of 
the signals is required prior to transmission. Because these 
Signals are unmodulated, only their relative position in any 
ten second time interval enables receivers to identify the 
Station Leing received. Jf primary importance in utilizing 
the OMEGA VLF signais are the ability to detect and measure 
the phase of a received signal, and the ability to predict 


the phase of the signai. Several factors have an effect on 
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tne propagation of VLF signals and conseguently the ability 
oI an airborne Omega system to receive and process the 
Sa IS ЕСЕ accutate positioming (Ref. 6]. These factors 
are described in the section pertaining to propagation 
Predierticn, 


2.  Receiver-Converter 


— s. ss h sss ue i s ы» s= mm ma m олы 


The AN/ARN-99 receiver was developed by tae 
Electronics Division O cuc Northrop Corporation min 
Houwrhrons, Calrftorznia. Іс was designed for operation with 
all eight transmitting stations and to provide a continuous 
update of aircraft Position: тһе receiver-converter 
installed in the P-3C consists of several modular assemblies 
as shown in Figure 6. The antenna, composed of two loop 
antennas fixed at right angles to each other and mounted at 
оаа е= to the aircraft center line, is controlled by the 
burst filter function in the computer through the antenna 
Switching matrix. The three heterodyne receivers have RF 
ama ІР sections loz processing of the input signals. The 
requirements for filtering the received signals in these 
sections are very stringent. The filtering in the RF 
section removes the RF image, including all frequencies in 
the агеа Of the harmonics of the local oscillator. The 
lari ow band filtering in the IF section provides for th 
ШОО ОПОТ ıimterferiıng carriers. In addition, limiter 
are used for controlling the dynamic signal levels, and a 
Dandpass filter is used to remove the harmonics created by 
the previous limiters. Ine resultant output of the RF and 
Ir filtering sections leaves the fundamental frequency as a 
Sinusoid. This output is fed into a sine and cosine 
Ccrrelator which utilizes a square reference signal to 
produce lan output or two direct current signals. One de 
EM ТЕ С proportional to tha sine of the burst phase, and 


the other dc signal is proportional to the cosine of the 
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burst phase. These dc signals, which represent the measured 
phase snift, enter the analog to digital converter and are 
transmitted to the conputer section through the 
receiver-converter/couputer interface box. The precision 
frequency generator provides for the reference signals and 
test signals necessary in the measuring of the phase 


difference of tne receiver signal [Refs. 7, 8]. 
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III. INIERNAL PROCESSING FUNCTIONS 


unap oe ee oe s D Í Ee e e He a амана ара a a E ae So St c ans 


The Omega navigation system in the P-3C aircraft is one 
of five navigation subsystens which comprise the  avionics 
navigation system. The remaining  subsystens, in their 
Eur LY OF Use, are a primary inertial, a secondary 
„ЫЫ а, a "doppler radar, and an air data computer. Tae 
dM econ and utilization of each subsystem is dependent on 
the control of tne Navigator/Comaunicator and оп the 
availability of the subsysten. A block diagram of the 


navigation subsystems is shown in Figure 7. 


A. OMEGA PROCESSING 


The internal processing functions OF the Omega 
navigation system are shown in Figure 8 (Ref. 4]. A more 
detailed description of the internal processing functions is 
presented in Figure 9 snowing the pertinent data transferred 


between the functions. 


Peers di iza tion 


Initialization of the Omega navigation system begins 
Ws the insertion of the date and time entries by the 
Navigator/Communicator and subsequent activation of the 
NM cuno top condition or marz aircraft position. The 
date and time entries are used within the propagation 
prediction module to calculate an estimate of the phase 


value derived from current aircraft position at the moment 
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FIGURE 8 - INTERNAL OMEGA PROCESSING FUNCTIONS 
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ence tjalızatıon, These phase values are used EO 
ze the tracking filters. The variance values within 
Re com national filter are also initialized at this tine, 
aba the position matrix R, utilized in the velocity and 
navigation processing function, is initialized based on 


wmr ent aircraft latitude and longitude. 


Hardware testing begins immediately after the 
termination of the initialization function, and also after 
саса restart. The main objective of the hardware testing 
module is ‚to determine the status от the Onega 
receiver/converter. Additionally, hardware testing providss 
a method of assuring proper operation of all Omega 
equipment. The testing involves a series of subtests wnich 
Malidates the operation or the input/output interface 
petween tne receiver/converter and the computer. The 
Subtests, which are described in greater detail in Ref. 4, 
are a communication subtest, a coherence status Subtest, an 
Omega output subtest, an Omega input subtest, a phase angle 
КО аста subtest,;, а phase counter subtest, and a RF/IF 
Ecc. An oscillator drift subtest, which is classified 
as a hardware test, is not implemented until single 
frequency processing ccmmences in the combinational filter 


routine, 


жн» қ) солы; MA GENS cee eee ош = ee 


БОЛСО ав оп Of the Omega navigation system to 
tae Omega transmission pattern is the most critical function 
executed by the navigation system. By aligning itself with 
the ten second interval transmission pattern, tae Omega 


navigation set processes the correct VLF signal bursts in 
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subsequent routines. Syachzennzaeıon 15 accomplished by 
utilizing the variable length transmission tine intervals of 
the varicus freguencies and a differential correlation 
technique. The first assumption made in the routine is 
that all eignt stations tnat are presently in operation are 
operating on all three frequencies. The routine, starting 
at some randca time, accepts from the receiver-cönverter a 
ten second interval of data composed of 100 sine and 100 
cosine values over contiguous 100 millisecond intervals. 
One hundred correlaticn coefficients are calcula:zed for the 

requency under consideration utilizing the difference 
between short sum and long sum intervals and adding these 
differences over the sight count burst pattern. Figure 10 
decicts the short sun and long sum intervals over a ten 


second interval pattern, and the limits of the short and 


long suas. The correlation coefficient formula is given 
below. 

8 1+JU 2 1+.JU 2 T+JU+2 2 JUE 2 
canes» > Xx (X)| + У 1 (K)| - У X (K)] * 2 Tal, 
J=1 | K=I+JL K=1+JL K=i+JL-2 K-ITJE-2 

(І- 1,100) JU-JL+1 JU-JL*5 


The values X(X) and Y(K) are the sine and cosine values 
respectively received from tne receiver-converter. Each of 
the confidence coefficients can be considered to bea 
neasurement cf how well the received ten second burst 
pattern matches the transmitted burst pattern. The higher 
Bnesererrieiene value, the greater confidence is associated 
ғо гПе index number I being the starting point of the 
transäitted burst pattern within + 50 milliseconds. The 
four largest confidence coefficients are saved in descending 
order of magnitude, along with their respective index number 
T. This index number is used in determining the starting 


Pommt O: the ter second burst pattern. Various conditions 
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Meee ug o the four highest confidence coefficients and 
index nunters are tested in order +o determine the success 
pu the failure of the synchronization attempt. ІЁ the 
current confidence level fails, the next sequential 
frequency is utilized to acquire the next ten second 
intervai cr data. The synchronization routine continues in 
Ес шійе until a süccessfiul synchronization occurs or until 


the Navigator/Communicator intervenes. 
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Шо ІІІ сег 


The burst filter routine conmences upon successful 
eomeletion of synchronization. This routine calculates the 
phase measurements and the variance of these phase 
measurements for each station and frequency received by the 
Omega havigation set. It also selects the optimal antenna 
pts and test data for the Omega receiver, maintains 
Ey menrG@gazZeation with the transmitted burst pattern, and 
ccaputes the phase measurement and its confidence 
measurement  Dased on the receiver phase difference and test 


data Iron tne receiver. The equations are: 


QJK - ARCTAN((XbJK)]|- OOK + O? 


(YbJK) 
Sr EM = QtX * (AtJ - QJK) 16 (02005) e е 


АСТ 2 (Ок = ок) 


where (0! is a phase correction (09 or 1809) added to the 
Be med signal to account or error in reception, OK is 
the phase angle ofíset for frequency K computed in the 
calibration measurement, and XbJK and YbJK are Durst 


measurement suns corrected for phantom error. 


To accomplish these calculations, each of the eight 
bursts and associated slots are divided into intervais of 
tage as shown in Figure il. The start burst and end burst 
intervals are periods of tiae where the sine and cosine data 
are not used. These intervais are considered transient 
emo ds Сг waiting periods. During the burst interval, the 


ШЕГІ сосіпес values received are used in the calculation 


th H 


S 
о the phase measurement and variance of the phase 


measurement. in the test interval, various measurements 
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are computed to improve the accuracy of the phase 
measurement, and for the issuing of antenna selection and 
test selection commands. reference 4, pages 86 to 101, 
describes in detail the computations involved in calculating 
the phase and the phase variance for each station and 


frequency used by the Omega navigation system. 
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Once the burst filter processing has acquired a 
phase measurement and an estimate of the variance of the 
Phase measurement for a specific station and frequency, 
G c values are used as inputs to only one of the 24 
m mn filters. Each filter is similar in the processing 
Ән he data, but each is different in its computational 
constants for the various frequencies and stations. A 
Same itied description of a tracking» filter is that it is 
utilized ЕС acquire a better estimate of the phase 
neasurement by averaging the outputs of the burst filter 
over successive ten second intervals. The tracking filter 
routine, using the measured phase from the burst filter and 
an updated estimate of phase based upon previous 
measurements and dead reckoning information, conbines the 
measured and estimated phase and utilizes a weighted average 
to compute a new phase estimate and phase variance estimate. 
Also computed within the tracking filter is the variance of 
the estimate of the error in the phase rate of change due to 
aircraft displacement in time. Once these three values are 
ccaputed for a particular station and frequency, a test is 
made to determine the accuracy of the phase measurement. Ir 
the variance of the pnase measurement is less than or equal 
to (0.06 fr radians)?, the pnase measurement is considered 
accurate enougn for computations within the combinational 
сег. After three successive measurements and successful 
Bests, a tlag is set indicating that the data is valia and 
available, and the combinational filter reads the outputs 
frcm the tracking iilter. Once the outputs have been read, 
the count of successive measurements is reset to zero, the 
variances of the phase measurements are set to the 
initialization values, and the data valid and available flag 


is cleared. If at any tine in the tracking filter the 
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variance cf the paase measurement exceeds (0.06 7 radians)?, 
the count of the successive measurements is reset to zero, 


and the data valid and available flag is cleared. 


6. Combinational Filter 


The input values to the combinational filter £ron 
tne tracking filter are the refined estimates of phase 
measurement, TRK, phase variance, o 2QTRK, and phase rate 
variance, o~ 2QERR. The inputs from the propagation 
prediction routine are the estimates of phase, (PROP, 
inproved for propagation effects and phase variance, 
ao 2PPROP. Within the combinational filter routine, these 
input values are statistically combined to provide the best 
estimates of systen position and velocity. The 
combinaticnal filter routine provides for the conversion 
from phase tc geodetic coordinates (latitude and longitude), 
and also determines the phase and frequency differences 
between the receiver oscillator and the transmitted OMEGA 
signal.  Cnce the Omega receiver oscillator is calibrated, 
the Omega navigation system operates when only two 
transmitting stations are accessible In addition, the 
combinational filter is used for lane determination through 


the technique of multiple state vectors. 


a. Description 


The combinational filter is a Kalman filter. 
Reference 4 describes the Kalman filter technique to 
filtering and prediction as a linear, recursive, minimum 
variance filter. R. G. Brown and L. L. Hagerman [Ref. 9] 
describe a Kalman filter as simply a means of estimating the 
various states of a random process from a set of discrete 


measurements having a known linear connection to these 
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States. It was further noted in Ref. 10, that a Kalman 
filter operates only on the system errors and not on total 
quantities such as position and velocity. The basic concepts 
involved are those of state, state transition, measurement, 
and optimal weighting [Ref. 11]. The states of the 
combinational filter are differentials of system variables 
ferme. CEEOL in position ($02 and $63), oscillator start time 
(t9), rate of change in oscillator start time ( EP and 
Serer іп velocity (dvc2 and бүсз) 1. The following chart 
describes the elements of the state vector X [6 x 1]. 


ELEMENT SYMSOLS MEANING 
1 642 Post cion error along the =R2 diracticn of 
the system position matrix R. 
2 603 Position error along the =-R3 directlcn of 
the system position matrix R. 
3 Ue Time difference between the transmitting 


station phase and the receiver oscillator 


phase. 
p Time rate of change of ts. 
óvc2 Error in the east component of velocity 


which is error along the R2 direction of 
the system position matrix R. 

6 бүсз Error in the north component of velocity 
САСА is error along the R3 direction of 


the system position matrix R. 


The state of the system is described by the 
solution of linear vector differential equations depicting 
system error growth. The Kalman filter method linearly 
combines the previous estimate of the state vector, X, with 
a measurement to approach a minimum variance estimation. 
This minimum variance estimate is then time updated until 
the following measurement. The measurement calculation is 


explained by the following equations. 
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Y = QTRK - QPROP 


where Y is the difference between the tracking filter phase 
measurement and the propagation prediction phase estimate of 


the phase measurement. 


Residual = Ү - M*X(t-1) 


where Residual (RES) is the difference between the 
calculated error in the measurement, Y, and the systen 
estimate of the error, H*X(t-1). The state vector X(t-1) is 
the previcus estimate of the state vector, and M[1x 6] is 


the measurement matrix described below. 


ELEMENT EQUATIONS MEANING 
3 
1 ЕЕ РЕ тж y R2J*SIJ EMER is the earths mean equa- 
J=1 torial radius = 20925741.47 LE 


K * SIN(DD) R2J and R3J are the elements of 
the R2 and R3 vectors of the 


3 position matrix R. 
2 EMER * » R3J*SIJ  SIJ's (12A,B,..,H) are the 
J= 1 elements of the station 


K * SIN(0D) direction cosine values. 
3 VL VL is the average propagation 
АК velocity, and AK (K=1,2,3) is 
the wave lengths of the 
frequencies 10.28,11.3, 13.6kHz 


3 
where фр = һассоз( 9 RIJ*SIJ] and the remaining elements of 
J=1 
M are zero. New, minimum variance, state vectors are 


produced by using the following equations. 
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änew = Xold + Б ж (RES) 
where b [6 x 1] is the optimal weighting vector matrix. 


This cptimal weighting vector, b, is established 
by means of time propagation of system error growth in 
combination with information Contained within the 
measurement. The weighting is a function of the covariance 
matrix, P (6 x 6], of the difference between system error 


State vector, X, and the average state vector Xavg. 
poc sse (X4- хата) * (X — Хауа)! 1 


where 2 is the expected value of the difference between the 
value of the difference between the value for the real state 
vector X, iZ all conditions чеге known, and the average 
state vector, Xavg. The following chart describes the 


elements of the covariance matrix P. 


PPE RED T SYTNEOLS MEANING 
11 o 202 The position variance for the positicn 
error, 602, along the -R2 direction 
P22 c 7403 The position variance for the position 


error, 663, along the -R3 direction 


P33 бг The variance of oscillator phase offset 

Puy 2 Thel vyatiance of Oscillator агы: rate 

P55 c ?VC2 The variance of the error in the east 
component of velocity óvc2 

F66 C ?VC3 The variance of the error in the north 


component of velocity ávc3 
All remaining elements of the covariance matrix are zero. 
The combinational filter effectively computes 


optimal estimates cf position and velocity along with the 


remaining elements of a state vector. The filter uses an 
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observation of the difference between the phase of the 
tracking filter and a value of the phase formulated on the 
time updated position of the combinational filter. Based on 
this measurement, the combinational filter derives estimates 
o£ positicn errors, velocity errors, oscillator start time 
ger OSClilator drift. 


Fecause the Omega navigation system functions 
with dy ral distinct velocity sensors (ie. inertial, 
doppler, or air data), a different node of operation exists 
for the combinational filter for each. sensor. Each method 
requires the filter to estimate, predict, and control a 
different set of system errors. Each set of system errors 
conforms to the dead reckoning velocity source used for tine 


updating. 


The combination of uncertainties in position 
БЕРОГ and receiver oscillator start time produces a 
vagueness of the number of wave lengths or lanes between the 
transmitting station and the receiver. This ambiguity occurs 
primarily when the system is first initialized. To find the 
solution to this lane ambiguity, several estimates of the 
E cEvector, X(i)'s, are calculated by the filter. These 
ШОН) "5 correspond to the different possible lanes or 
integral values of phase measurement. As measurements are 
acquired, corresponding to the three frequencies, the 
Uuscertalnity in oscillator Start time, tg, is diminished. 
Measurements from different stations tend to minimize the 
position uncertainity. Ihe general effect is a reduction in 
the number of state vector estimates, X(i)'s, that can be 
considered logical estimates of the most accurate state 
vector until only one credible estimate remains. The 
Gemperion for reasonableness is based on variance 


considerations. 
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D. Operation 


There are three distinct operations identified 


Bach the combinational filter routine. 


© En2tıalization 
2) Tine update 


3) Measurement update 


The initialization operation assigns the initial 
values of variance to the covariance matrix P, sets the 
elements cf the initial state vector to zero, and sets the 


system driving noises as a function of the velocity source. 


In tne time update operation, both the state 
vector X, and the covariance matrix P, are predicted from 
the last update. The following equations are utilized in 


the time update operation. 


Xnew - Ónew * Xold 


= T 
Pnew = Önew * Pold * Qnew + Nnew 


The transition matrix ф [6 х 6] mathematically expresses the 
Beesagaeıon of errors across the time interval At since the 
last update. The elements of the transition matrix are 


shown in the following chart. 


1 0 0 At 0 
0 1 0 0 At 
0 0 1 At 0 0 
0 0 0 0 0 
0 0 0 1-Bi*At 0 
0 0 0 0 1-gi*At 


where Bi (i=T,D, or à) signifies the 1пуегзе correlation 
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time of бүс2 апа $vc3 for the velocity source mode utilized. 
The transition matrix f furnishes the propagation of the 
predicted elements of the state vector X and the covariance 
matrix P across the time interval At. The variance effects 
across this time interval are propagated by the additive 
diagonal псісе matrix N [6 x 6]. The following chart depicts 


the elements of the noise matrix N. 


0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 Biro2ixAt 0 

0 0 T9 0 0 Qi*oc?i*At 
Biers Bi is the inverse correlation tine described 


previously, and c 21 is the variance of the velocity errors 
(óvc2 and óvc3) and is dependent on velocity source utilized 


Os D. or M. 


During the time update operation, the multiple 
state vectors, Х(1)!5, which exist until the correct lane 
has been determined must also be time updated. The time 
update equation for the state vector X must be sequenced 


through ali state vectors. 
X(i)new = Qnew * X(i) old fOr ad a 


The covariance matrix time update equation is computed once, 
because there is only one covariance matrix regardless of 


the number o£ state vectors. 


After completion of the time update operation 
and prior tc the measurement update operation, the 
measurement residual, (RES), measurement matrix, M, and the 


measurement confidence scalar, C, are calculated for each 
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tracking rider input. The measurement residual and 
measurement matrix caiculations were previously explained in 
the e e on describing the Combinational filter. The 
measurement confidence scalar, C, is a measurement of noise 
computeđ by the following equation. 


С = с 2$PROP + c 2фТЕК 


where o 2(PROP and o- 2OTRK are the variance of 
the propagaticn prediction phase estimate and the tracking 
filter phase measurement. 


The measurement update operation involves the 
generaticn of a linear, unbiased weighting vector b, and 
computation cf a new, minimum variance, state vector X by 


the formula listed below. 
Xnew - Xold * b * (RES) 


TO generate the optimum weighting vector, b, three 


expressicns are required. 


1) the predicted measurement variance - мж ржи! 
2) the measurement confidence scalar - C 
3) the divergence control factor - € (epsilon) 


The weighting vector is calculated from the following 


equations. 


мжржи! ЕКЕ 


tO 
II 


(Dem кезі) 2 


o 
II 


where Q is defined to bea scalar quantity depicting the 
IEEE riasnce, and € is a factor to prevent divergence of 
the estimaticn process. Also updated in the measurement 


update process is the covariance matrix based on the 
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following formula. 


Pnew - Pold - b * M * Pold + € x м! ж gj 


The retention of new state vectors is predicated on the 
testing cf the position elements in the state vector, and 


the variance values in the covariance matrix. 


Jee eePpagation Prediction 


In crder to effectively resolve the aircraft 
position frcn the phase information received, the phase 
velocity cf the wave along the propagated path and tne 
length of the path must be known. Several factors or 
effects create imperfections in the wave guide through which 
the VLF signals propagate, The propagation prediction 
function establishes the best estimate of phase and phase 
variance taking into account the factors listed below which 


effect the phase velocity of the propagated signals. 


a. Diurnal Effects 


The positionai changes of the sun over the earth 
adjust the size and the shape of the ionosphere. The change 
of the wave guide through which the VLF signals propagate 
causes a change in the velocities of the phase signal which 
results in a difference in phase angle received at the Omega 


receiver location. 
b. Ground Conductivity 


This factor takes into account the effect o£ 
various mediuns over which VLF signals travel. Water, which 
EN pperfect conductor of VLF signals, does not greatly 
affect the signals. Land, which is a less perfect conductor 


of VLF signals, does affect the signals to a greater extent. 
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Conductivity patterns have been measured and are predictable 
Ier. 6]. 


C. Earth's Magnetic Field 


This factor affects both the attenuation rate 
End спе velocity of tke signal. The phase shift on 
refracticn of the VLF signal off of the ionosphere depends 
спе іпгегасізлоп ос <Һе signal with the electrons, and 
this depends on the orientation of the magnetic field with 
respect to the direction of propagation and on the magnitude 
cri ne field [Ref. 12]. 


d. Latitude and Spheroidal Effects 


This factor takes into account the nonspherical 
Shape of the earth and the adjusted path of travel from the 


transmitting staticn to the Omega receiver. 
e.  Folar Cap Absorption (PCA) 


VIP signals which pass over the Polar Сар 
experience an abnormal rate of absorption during both day 
and night. PCA produces large changes in VLF signal patterns 


which can last for several days. 


- 


|. f. Sudden Ionospheric Disturbance (SID) 


Sclar flares emanating on the sun's’ surface 
increase the ionization within the atmosphere. These SID's 
cause large changes in the phase of VLF signals passing 


through the area of activity. 


Of the six factors which affect the propagation of 
ҮЕ Signals, the first four are predictable and can be 
compensated for algorithmically, while the last twc are 
unpredictable and cannot be compensated for within an Omega 


navigation system. 


43 





The velocity and navigation processing function is 
divided into two routines in the Omega navigation systen. 
The velccity processing routine, utilizing the current 
navigaticnal source (ie. inertial, doppler, or air data), 
calculates the velocities along the system axes frcm the 
available scurce velocities. These system axes velocities, 
conbined with the state vector corrections from the 
combinational filter, are supplied to the navigation routine 
in order to update the Onega systen position matrix and the 
Omega aircraft latitude and longitude position. 


IE LOCITY PROCESSING 


ИП ше velocity processing routine, the aircraft current 
true airspeed velocity (VTAS) and heading (ACHDG), the 
current wind direction (WD) and velocity (VW), and the 
dopplers current velocity components along heading (VDA) and 
across heading (VDC) are utilized to compute the velocities 
in the North/South and East/West directicns from the 


following equations. 


VEN = VDA * COS (ACHDG) VDC * SIN(ACHDG) 
VDE = VDC * COS(ACHDG) + VDA * SIN(ACHDG) 
VAN = VTAS* COS (ACHDG) + VW * COS(WD) 
VAE = VTAS* SIN(ACHDG) + VW * SIN(WD) 


where VDN and VDE signify doppler velocities North/South and 


East/West, and VAN and VAE signify air data velocities 
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North/South and East/West respectively. In constructing the 
data flow graph of these equations, it was self-evident that 
some variables could be calculated, stored in nenory, and 
retrieved from memory when they were required in later 
calculations. Figure 12 depicts the data flow graph for the 
North/South and  East/West computations. Implementing PL/M 
code directly frcm this data flow graph was straight 
forward. 


On completion o£ the Nortn/South and East/West 
velocities computations, a test is made in the velocity 
processing functicn to determine if the navigation mode has 
changed since the last iteration or the velocity function. 
If an upgrade of the navigation mode has occurred (ie. 
doppler tO inertial, or air data to doppler or inertial), 
then the summations (DELVC2 and DELVC3) of the corrections 
to the system velocity components from the combinational 
filter are set equal to zero, and new velocities along the 
system axes are calculated. Figure Gepicrs пе 
navigaticn mcde test and applicable equations for computing 


the system axes velocities. 


ІП Ter го Save execution time am the caiculation of 


the formulas and to simplify the PL/M code, the functicns 


(VÍE * COS(SHDGA) + ViN * SIN(SHDGA)) 
(VİN * COS(SHDGA) - ViE * SIN(SHDGA)) 


were calculated prior to the navigation mode test. Figure 14 
deriets the data flow graph for these two functions. After 
successful computation of the above functions and the 
navigaticn mode test, the corrected velocities (VC2 and VC3) 
alcng the system axes are calculated and the resultant 


‘values used in the navigation routine. 
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VDA :k COS(ACHDG) 
> 
VDA = VDN 


n * SIN(ACHDG) 
ACHDG |COS LA 
SIN 


| УОС ж COS(ACHDG) 
E '* SIN(ACHDG) 


Ке * COS(ACHDG) 


VTAS VAN 
-- 


WD cos 
SIN 
VTAS * SIN(ACHDG) 
ыт 
VW VAE 
EIC * SIN(WD) 


FIGURE I2 - NORTH/SOUTH/EAST/WEST VELOCITY DATA FLOW GRAPH 
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В. NAVIGATICN PROCZSSING 


In the navigation processing routine, the Omega system 
POSICIONES RR [3 X 3] is initialized, the angular 
rotations  (DELT2 and DELT3) about the system axes are 
calculated using the corrected velocities supplied by the 
velccity processing routine and the state vector corrections 
seemwmthe cefbinational filter, tne system position matrix is 
updaced, and a new Omega aircraft latitude (OLATAC) and 
longituds (OLONAC) are calculated. The flowchart of the 


ПӘМЕСасісп routine is shown in Figure 15. 


The initialization of the Omega system position matrix R 
occurs during the initialization function when the 
Navigator/Communicator inserts the aircraft initial latitude 
(OLATIN) and longitude (OLONIN). These values, with an 
initial ccndition that the system heading angle (SHDGA) is 


zero, are input to the equations listed below. 


R11 = SIN(OLATIN) 
R12 = COS(OLATIN) * COS(OLONIN) 
R13 = COS(OLATIN) * SIN(OLONIN) 
P21 = COS(OLATIN) * SIN(SHDGA) = 0 
R22 = -SIN(OLONIN) * COS(SHDGA) - 
SIN(OLATIN) * COS(OLONIN) * SIN (SHDGA) 
= -SIN(CICNIN) 
R23 = COS(OLONIN) * COS(SHDGA) - 
SIN(CLATIN) * SIN(OLONIN) * SIN (SHDGA) 
= COS(CICNIN) 
R31 = COS(OLATIN) * COS(SHDGA) = COS (OLATIN) 
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R32 


SIN(OLONIN) * SIN(SHDGA) - 
SIN(OLATIN) * COS(OLONIN) * COS (SHDGA) 
-SIN(CLATIN) * COS (OLONIN) 


233 = -COS(CLONIN) ж SIN(SHDGA) - 
SIN(CLATIN) * SIN(OLONIN) * COS (SEDGA) 
-SIN (OLATIN) SIN(OLONIN) 


+ 


These nine equations represent nine direction cosine values. 
These values represent the orientation of each of the three 
Euteuuwxes (RiwR2, and 83) in the earth's fixed coordinate 
frame. E MEUCOUuStrücting the data flow огарие lor the 
wnk alization equations, the concept OL Mininizing 
execution time prevailed, and as variabies were calculated, 
they were stored in memory and retrieved as they were 
Secured 1n later calculations. Figure 16 portrays tae data 


ME rapa icr thé initialization of the position matrix R. 
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FIGURE I6 - INITIALIZATION OF POSITION MATRIX R 
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After completing the Initialization phase, the 
Beeren Zunetion calculates the angular rotations (DELT2 
and DELT3) about the system axes R2 and R3. The following 


L Ons are used in calculating the angular rotations. 


БЕО Е ЕТЕ EEC = (2 = R312 - 8112) ] ж УСЗ + DINAY > 
EMER 
2 * EEC * R21 * R31 * VC2 * DTNAV + 663 
EMER 
В ИІЗ - [1 + EFC &-(2 * 2212 - к11г) | * YC2 * DTNAV + 
^— BHER 
2 Ж БЕС * R21 * 331 * VC3 * DTNAV - 602 


ENER 


As in previous equations, converting to data flow graphs was 
accomplished with the idea of minimizing execution time. 
Figure 17 is the data flow graph constructed from the above 
emuations. Once the angular rotations are calculated, the 


position matrix is updated based on the formuia 


[Alnew 2» (R]old + [$] * [RJold 


where [0] is the rocation update matrix given below. 


- (DELI22 + DELT32) DELT3 -DELT2 
2 
-DELT3 - (DELT3) DELT2 * DELT3 
2 2 
DELT2 DELT2 * DELT3 - (DELT2) 
2 2 


simplificaticn oí the position matrix update formula 15 


accomplished utilizing the following equations. 
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(j) = ij + R25 * (DELT3 7 2) - R3j * (DELT2RV 2) 


B = [ (DELT2)2 + (DELT32) ] 72 


The position matrix update equations are then written in the 


following format. 


R1j = R15 - R13 * B + R2j ж DELT3 - R3j * DELT2 
R2j = R2j - DELT3 * A(3) 
BOE SER3j * DBLT2 *. A (3) 


Inplementaticn of these equations in a data flow graph is 


HSpaoted in Figure 18. 


Once tte position matrix R is updated, the last 
equations to calculate are for updating the system heading 
angle and the aircraft new latitude and longitude. The 


equacions listed below update these three system variables. 


SHDGA = ARCTAN (R21 / R31) 
OLATAC = ARCTAN [ (R11 * COS(SHDGA)) / R31] 


ARCPAN (R13 / R12) 


OLONAC 


The data flow graph presented in Figure 19 depicts the 


update of these system parameters. 
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МЕСІ, che data Zlow graphs of the velocity and navigation 
function, four PL/M floating point subroutines were written. 
The subrcutines consist of an initialization (INITIALIZESR) 
routine for initializing the position matrix R, a velocity 
processing (VELPROC) routine for computing the velocities 
along the system axes, a navigation processing (NAVPROC) 
routine for calculating the angular rotations about the 
system axes, and a new position (NEWPOS) routine for 
updating the position matrix R, aircraft system heading 
angle, and aircraft latitude and longitude. Based on the 
arithmetic and/or trigonometric functions used in the data 
flow graphs and an approximation of the time required to 
execute each function, a preliminary estimate of the time 
reguired to execute each subroutine was acquired. Utilizing 
INTERP/80, which is the PORTRAN IV software simulation of 
the INTEL 8080 CPU and is available on the IBM 360/67, the 
actual execution time of each program was acquired. The 
following chart shows the number of arithmetic functions 
used in each subroutine, the estimated and actual execution 
times in milliseconds of the subroutines, and the storage 


utilized for each subroutine. 


ADD SUB MULT DIV COSSIN AT AN 
(1.0-1.2) (1.5-2.0) (55.0-66.0)  (35.0-40.0) 
INITIALIZESR 0 0 4 0 2 0 
VELPROC 6 6 12 0 3 0 
NAVPROC 3 3 11 0 0 0 
NEWPOS 10 9 18 3 1 3 


ESTIMATED “ACTUAL T STORAGE 


TINE DENE (BYTES) 
INITIALIZESR 123.9 113.4 375 
VELPROC 206.7 248.4 830 
NAVPXOC Dan Zo 384 
NEWPOS ЛЕТ а 842 
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Fron the last three subroutines, it is seen that the total 
execution time of 562.5 milliseconds is between two and 
one-half and three times as great as the 200 millisecond 
time interval currently implemented in the P-3C. However, 
because the Omega navigation display for Omega latitude and 
longitude is updated every second, the microcomputer 
subroutines are capable of providing an update in sufficient 
time to be displayed on the Omega tableau. 


In the simulated test case used to check the validity of 
the subroutines, a navigation cycle time (DTNAV) of 1 second 
was used in the computations. This 1 second time interval 
allowed an aircraft traveling at 300 knots to be updated at 
approximately every 506 feet. The velocity processing 
routine was able to calculate the velocities along the 
ENStem faxes, using the PL/M floating point function, to four 
digit accuracy (ie. from 000.0 feet per second to 999.9 feet 
per second). These velocities encompass the complete range 


of the P-3C airspeed capabilities. 


The calculations performed in the remaining subroutines 
presented problems when using the PL/M floating point 
functions.  Eecause the functions utilize a 16 bit mantissa, 
the exponent values have to be within a range of 216 when 
adding or subtracting two numbers. Because the difference 
between the calculated angular rotations about the system 
axes and the direction cosine values in the position matrix 
is not within this range, several angular rotations 
calculated cver time had to be added together before 
updating of the position matrix R could be accomplished. 
This is the major difficulty with the subroutines 


implemented with the present PL/M floating point routines. 
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Based on the subroutines implemented for the velocity 
and navigation function, a microcomputer is capable of 
calculating the required system parameters in sufficient 
time and with sufficient accuracy to provide reliable 
position fixing information. Even with the microcomputers 
primary disadvantage of slower execution time, the position 
fixing information could be displayed at one second 
intervals provided a floating point package which operated 
with a mantissa of 24 bits could be acquired. Further 
research of the remaining internal functions and the overall 
effect they have on the velocity and navigation function is 
encouraged. Particular attention should be focused cn the 
State vector elements applied to tne velocity and navigation 
equations, and their relevance to overall system operation 


апа update of positicn information. 
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КОБЕ 77 ACCRESS, ZE BYTE: 
/ = 
THE ABCVE VARIABLES ARE GLOBAL 
TO ТРЕ FLOATING PCINT SUBRGUTINES 
ERRCF: FROCEDURE(I); 
/* 
PRCCECURE TO PROVIDE ERROR MESSAGE 
C OLLCWED RY TERMINATION CE PROGRAM EXECUTION 
DECLARE I BYTE: 
CALL CRLF; 
ie CASE wiles 
/> 
MESSAGE FOR UNDERFLOW 
CALL PRINT(.'EXECUTICN TERMINATEC DUE TO UNCERFLCM$'); 
/* 
MESSAGE FOR OVERFLOW 
* 
CALL PRINT(.'EXECUTION TERMINATED DUE TO CVERFLOWS'); 
/ 
MESSAGE FOR ERROR IN LCG-ARGUMENT 
CALL PRINT(.'ARGUMENT FCR LOG NEGATIVE QR ZERCS'); 
/* 
MESSAGE FOR ERROR IN OIVISION 
CALL PRINT(.'ATTEMPTED DIVISION BY ZERO$S'); 
/5 
MESSAGE FOR INPUT ERROR 
/ 
E OMISPRINTU.UINPUT ERRORS); 
/% 
MESSAGE FOR EXP(X) OVERFLOW 
/ 
CALL PRINT(.'OVERFLOW, ARGUMENT FCR EXP(X) TCC BIG£'); 
END; 
FALT; 
RETURN; 
END ERRCR; 
/ж Д 
ELCATING POINT ACD ROUTINE 
= 
ADO: PROCEDURE (ХА,ҮА,РА); 
CECLARE XÀ ALCRESS, X BASED XA BYTE, 
YA ACLRESS, Y BASED YA BYTE, 
PA ADDRESS, P BASED PA BYTE, 
UDBEEST. SIGN, XEX,YEX) BYTE; 
(XX,YY) ADDRESS; 
CECLARE FOSITIV LITERALLY '< 80H', 
REGATIV LITERALLY '=80H'; 
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PROCEDURES 





DECLARE 


H 


LL CCSSIN (.ACHDG, .В8, .А); 


Ж Rn 


СС5(МІМО ПІКЕСТІСМ), С - СІМ(МІМО DIRECTION) 
EE CESSIN (205 Ds, C); 


SINCACHDG) * VDC 
LL NULT (.A5, .VDC, .А1); 


z CCS(ACHDG) * VDAÀ 
LL MULT (.B, .УПА, .В1); 


PN AB Dor D 


= SIN(ACHDG) * VDA 
LL MULT (.A, .VDA, .A21; 


OPN ODN ODN С ЗЕМ OO + 


І» “У №) 4 


CESPACHDGI) + VIC 
LL MULT UND, VDC; БЕЛІ 


IX OIDO 
O РММ 


ET 


= COS(ACHDG) ж VDA - SIN(ACHDG) * VDC 
LL SUB "Bl, SAT, “УОМ); 


= СОЅ(АСНОС) * VCC + SIN(ACHDG) * VDA 
L ADU (.B2, .А2, -VDE); 


~ m 


CCS (ACEDG) + VIAS 
CE MULT ЕВ. SUMAS 292): 


SIN(ACHDG) * VTAS 
LU MULT GA; ЫМ До» „АЗ3}ў 


CGS(wl) * VW 
LL MULT р, «YA; «DI 


- SIN(WE) * VW 
LL NULT (.C, .VW, .C110; 


- CCS(ACHDC) * VTAS + COS(WD) ж VW 
L ADD (.B3, «Dl, VAN); 


Okkes СОнсс- СС. Оғы. От. Ог<7- MH 
tM AR Dm Deae 4 Duo J a DOOR D 


[x 


85 


CCS(AIRCRAFT HEALINC); A = SIN(AIRCRAFT HEADING) 





/* 
VAE 
ж 


Y = 


) SIN(ACHDG) 
CALL AED (.A3, 


whe 
Ts 


СІ, 


; CCS(SHEGAI, A = 
BERSESCESSIN I .SHDGA, 


VTAS * SIN(WD) 


VAE); 


S IN( 


SHDGA) 


* VW 


COMPONENTS =/ 


COMPONENTS */ 


CCMEONENTS */ 


„8, +A); 

CC CASE NAVMODE; 

/* CASE О = MNERT PAL 
VELN(0) = VIN(0): 
VELN(1) = VIN(1); 
VELN(2) = VIN(2)3 
VELE(O) = VIZ(0)5 
VELE(L) = VIE(1); 
УЕГЕ(2) - УІЕ(2); 
END; - 

/* CASE 1 = DOP EI ER 
VELN(O) = VDN(C); 
VELN(L) = VDN(1)3 
VELN(2) = VDN(2)3 
VELE(O) = VDE(O); 
VELE(1) = VDE(1); 

ЧЕГЕ (2) = ҮПЕ(2); 
END; 

/* CASE 2 - AIR DATA 
VELN(2) s» VAN(GO); 
МЕСМ(1) = МАМ( 1): 
VELN(2) = VAN(2); 
VELE(3)) = VAZ(O); 
VELE(1) = УАЁ(1); 
VELS(2) = VAE(2); 

END; 

ENC; 

/* 

81 ЕМЕШ *-COSCSPOGA) 

/* 

Al = VELN * SIN(SHOGA) 

ж / 

Came gtr (.VELN, -A7 Al) 

/ 3 

22 = VELN * COS(SFOGA) 

* / 

(АО DULT (VELN+ .B, .в2); 

[> 

А2 = VELE * SIN(SFDGA) 

х 

MOT УЕГЕ, 245.42): 

/* 

A2 = VELE * COS(SFOGA) 

= 7 

CAEL AED (BI, ЖАП» САЗ): 

/* 

RS = VELN * COS(SHOGA) 

* 7 

CALL SUB (2.22, 55225 „В3); 
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+ VELN * 
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- VELE * SIN(SHDGA) 


ОР МАМ. MODE 
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